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We investigate to what extent a description of 1 Be as a three-body system made of an inert 
10 Be-core and two neutrons is able to reproduce the experimental 12 Be data. Three-body wave 
functions are obtained with the hyperspherical adiabatic expansion method. We study the discrete 
spectrum of 12 Be, the structure of the different states, the predominant transition strengths, and 
the continuum energy spectrum after high energy fragmentation on a light target. Two + , one 2 + , 
one 1~ and one 0~ bound states are found where the first four are known experimentally whereas 
the 0~ is predicted as an isomeric state. An effective neutron charge, reproducing the measured 
B(E1) transition and the charge rms radius in 11 Be, leads to a computed B(E1) transition strength 
for 12 Be in agreement with the experimental value. For the EO and E2 transitions the contributions 
from core excitations could be more significant. The experimental 10 Be-neutron continuum energy 
spectrum is also well reproduced except in the energy region corresponding to the 3/2~ resonance 
in n Be where core excitations contribute. 

PACS numbers: 21.45.-v, 21.60.Gx, 31.15.xj, 27.20.+n 
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I. INTRODUCTION 

The second lightest bound nucleus in the N = 8 iso- 
tonic chain is 12 Be, placed in the nuclear chart just above 
the widely investigated borromean n Li nucleus. These 
two nuclei are essential to understand the breakdown of 
the N = 8 shell closure when the dripline is approached. 
The parity inversion in n Be, already known in the early 
70's [J], is a clear indication of this fact. The particle 
unstable nucleus 10 Li should in principle have the same 
neutron configuration as 11 Be. Several theoretical works 
also predicted the existence of an intruder low-lying s- 
wave state @, S S|- The available experimental data 
concerning the ground state properties of 10 Li are how- 
ever controversial, although most of them point towards 
the existence of such low- lying virtual s-state [a, fi M B • 
This result has been confirmed in the recent work [9j . 

These properties of n Be and 10 Li clearly suggest that 
the p — sd shell gap is reduced when approaching the 
neutron dripline, leading to a structure of 12 Be and n Li 
with a large contribution from sd configurations. This 
has been confirmed theoretically and experimentally for 
both 12 Be [IS EH IH El and n Li 0, Q El] . 

The n Li properties have been successfully described 
by use of three-body models that freeze the degrees of 
freedom of the 9 Li core [lg, [l7|- It is then tempting to 
follow a similar procedure to investigate 12 Be. In fact, 
although 12 Be is not borromean, 11 Be is considered to be 
the prototype of one-neutron halo nuclei [18l . Il9l . |2(J ] , and 
therefore a description of 12 Be as a 10 Be core surrounded 
by two neutrons appears as a good first approach. 

An important advantage of 12 Be compared to n Li is 
that the properties of n Be are are much better known 
than the ones of 10 Li, and therefore the uncertainties aris- 
ing from the core-neutron interaction should in principle 



be smaller. However, while in n Li the 9 Li core is spher- 
ical, the 10 Be core in 12 Be is deformed, and as such one 
of the essential reasons for the shell closure breaking in 
12 Be. The ground state in n Be contains an important 
contribution from core excited configurations [21]. Dif- 
ferent theoretical calculations have estimated this contri- 
bution, and the results range from 40% in [22| to 10% in 
[Hij, passing through 20% in (HHf!. Recent experi- 
mental data [27, 28] are consistent with a 16% admixture 
of core excitation in the n Be ground state wave function. 

As a consequence of this, a description of 12 Be as an 
inert core surrounded by two neutrons is quite question- 
able, and the role played by the core excitations is an 
important issue to be clarified. In [291 ] a hyperspheri- 
cal expansion of the three-body wave function was used 
to obtain the ground state of 12 Be including core exci- 
tations. It was found that simultaneous fitting of the 
experimental ground state energy and the experimental 
longitudinal momentum distribution of 10 Be after high 
energy fragmentation of 12 Be, required a strong core ex- 
cited component in the wave function (sa 42%). 

However, after publication of this work additional ex- 
perimental information about the spectrum of 12 Be be- 
came available. A bound 2 + state and a bound 1~ state 
were found with excitation energies of 2.10 MeV [30ll3l| 
and 2.68 MeV [32J, respectively. Also the existence of 
a second + bound state with excitation energy of 2.24 
MeV was already envisaged [33| (later on confirmed [34j]). 
Also, the ( 12 Be, 11 Be2J one neutron removal measure- 
ments at the NSCL [12J allowed a direct estimate of the 
( 11 Bc(j 7r )| 12 Be(gs)) spectroscopic factors. All this new 
data lead the authors of [29| to review their calculations 
[35| . They found that a reasonable simultaneous match- 
ing of the data required a significant reduction of the core 
deformation, and hence a smaller contribution from core 



excitation (« 20%). 

Very recently new experimental data on 12 Be have 
been provided [13[. Specially interesting is the contin- 
uum relative energy spectrum for 10 Be+n after high en- 
ergy breakup of 12 Be on a carbon target. This invariant 
mass spectrum is known to be very sensitive to the fi- 
nal state interaction [36l |37|. and therefore, in our case, 
to the properties of the unbound 11 Be states. A re- 
liable calculation of the spectrum requires inclusion of 
the 10 Be-neutron continuum states together with core- 
neutron resonances. This energy spectrum is then a very 
useful observable to investigate the role played by the 
10 Be-neutron interaction, and to constrain the remain- 
ing uncertainties in the structure of 12 Bc. 

In previous 12 Be three-body calculations ( 10 Be+n+n) 
[J, [23, [35| the wave functions were obtained using the hy- 
pcrharmonic expansion. The employed neutron-core po- 
tentials were chosen to reproduce the bound n Be states 
and perhaps the first resonance at 1.78 MeV (excitation 
energy) but higher unbound states were ignored. This 
method is not the most efficient for a non-borromean 
system like 12 Be, where two of the two-body subsystems 
have bound states. This is because an infinite hyper- 
harmonic basis is in principle needed to reproduce the 
correct two-body asymptotics [38|. The convergence of 
the hyperharmonic expansion is slow, and in practice the 
energies and rms radii are obtained after extrapolation 
of the numerical results [J, [23, [35| where the basis is 
progressively increased up to a maximum value of the 
hypermomentum K max = 20 for the ground state, and 
K max = 12 for the excited states. 

In the present context the hyperspheric adiabatic ex- 
pansion method [391 ] is more appropriate. This method 
solves the Faddeev equations in coordinate space, treat- 
ing symmetrically the three two-body interactions such 
that each of them only appears in its natural coordi- 
nates. This makes the method specially suitable for non- 
borromean systems like 12 Be, where more than one two- 
body subsystem has bound states. Also the method per- 
mits the use of much larger values of the hypermomen- 
tum quantum number which guaranties convergence of 
the results. 

The main goal of the present work is to assess to what 
extent a three-body model with an inert 10 Be-core is 
able to reproduce the existing rather large amount of 
both old and new experimental data concerning 12 Be, i.e. 
the 12 Be spectrum, the El and E2 transition strengths 
[HI, [H, H|, the MCEO) [11 and the measured invari- 
ant mass spectrum [13J. The two-body potentials to be 
used in the calculations are fitted to reproduce the avail- 
able two-body experimental data, in particular, the Be 
data. In this sense, although the 10 Be core is considered 
an inert particle with spin and parity + , the employed 
two-body interactions phenomenologically account for all 
effects of core excitation appearing in the corresponding 
channel. In particular it is interesting to know if the 
weakly bound excited states can be understood as halo 
states and perhaps can be better described in a three- 



body model than the relatively well bound ground state. 
Relations between various quantities can be tested and 
new properties predicted. For this the hyperspheric adi- 
abatic expansion method is well suited. By comparing 
computed results and available data we can establish how 
large the contributions must be from inclusion of core ex- 
citations. 

The paper is organized as follows. In section |TT] we 
very briefly describe the basis of the three-body method. 
The different two-body interactions used in the calcula- 
tions are detailed in section IIII1 The results are shown 
in sections IIV1 fVl and IVI1 where we discuss the spec- 
trum and structure of 12 Be, the electromagnetic tran- 
sition strengths, and the invariant mass spectrum after 
high-energy breakup, respectively. We finish in section 
IVIII with a summary and the conclusions. In the ap- 
pendix some remarks about the El and E2 operators 
are given. 



II. THEORETICAL FORMULATION 

We assume 12 Be can be described as a three-body sys- 
tem made by a 10 Be core and two neutrons. The wave 
functions for the different bound states are obtained with 
the hyperspherical adiabatic expansion method. A de- 
tailed description of the method can be found in [3!| • 

This method solves the Faddeev equations in coordi- 
nate space. The wave functions are computed as a sum of 
three Faddeev components rp^\xi, yj (-£=1,2,3), each of 
them expressed in one of the three possible sets of Jacobi 
coordinates {a^,?/,}. Each component is then expanded 

in terms of a complete set of angular functions {<fi„ } 
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where p — \j x 2 + y 2 , a% — arctan(xi/yi), Vt Xi , and Cl Vi 
are the angles defining the directions of {a:.; and y^\. 
Writing the Faddeev equations in terms of these coor- 
dinates, they can be separated into angular and radial 
parts: 
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where Vjk is the two-body interaction between particles 
j and k, A 2 is the hyperangular operator [391 ] and m is 
the normalization mass. In Eq.Q E is the three-body 
energy, and the coupling functions P nn > and Q n n' are 
given for instance in 39]. The potential V^(p) is used for 
fine tuning to take into account all those effects that go 
beyond the two-body interactions. In the present cases it 



is rather small and unless the opposite is explicitly said, 
this three-body potential is taken equal to zero. 

It is important to note that the angular functions used 
in the expansion |T]) are precisely the eigcnfunctions of 
the angular part of the Faddeev equations. Each of them 
is in practice obtained by expansion in terms of the hy- 
perspherical harmonics. Obviously this infinite expan- 
sion has to be cut off at some point, maintaining only 
the contributing components. 

The eigenvalues \ n {p) hi Eq.0 enter in the radial 
equations © as a basic ingredient in the effective radial 
potentials. Accurate calculation of the A-eigenvalues re- 
quires, for each particular component, a sufficiently large 
number of hypcrspherical harmonics. In other words, the 
maximum value of the hypermomentum (K max ) for each 
component must be large enough to assure convergence of 
the A-functions in the region of p- values where the f n (p) 
wave functions are not negligible. 

Finally, the last convergence to take into account is 
the one corresponding to the expansion in Eq.|T|). Typ- 
ically, for bound states, this expansion converges rather 
fast, and usually three or four adiabatic terms are already 
sufficient. 



III. TWO-BODY INTERACTIONS 



It is known that for weakly bound systems, like for 
instance 6 He or 11 Li, the short distance behaviour of the 
two-body potentials is relatively unimportant as long as 
they reproduce the low-energy scattering data. Then the 
essential three-body properties can be described [la ]. 



A. Neutron-neutron potential 

For the neutron-neutron interaction we use a simple 
potential reproducing the experimental s- and p-wave 
nucleon-nucleon scattering lengths and effective ranges. 
It contains central, spin-orbit (£ • s), tensor (S'12) and 
spin-spin (si • s 2 ) interactions, and is explicitly given as 
E3 

^(rHST.OSe^/ 1 - 31 ) 2 ) 

+ (49.40e(-(''/ 1 - 31 ) 2 ) + 29.53e(-(''/ 1 - 84 > 2 >) Sl ■ s 2 
+7.16eS- (r / 2A3 ^S 12 , (4) 



where £ is the relative orbital angular momentum be- 
tween the two neutrons, and s = Si + s 2 is the total spin. 
The strengths are in MeV and the ranges in fm. We will 
refer to this potential as gaussian neutron-neutron po- 
tential. 

To test the role played by the short-distance proper- 
ties of the neutron-neutron interaction, for some specific 
cases, we are also using the more sophisticated vs ver- 
sion of the nucleon-nucleon Argonne potential [42| . This 



TABLE I: Strengths (in MeV) and ranges (in fm) of 
the central and spin-orbit gaussian potentials (V c (r) = 
S W e -(^ } )\ V W {r) = sWe-cMlV ) entering in Eq.© 
for the four interactions used in the calculations. For d-waves 
the radial potentials are made as the sum of two gaussians, 
whose strengths and ranges are given by the corresponding 
two rows in the table. In Potential IV the acronym "P.E.P". 
refers to the Phase Equivalent Potential used in this case for 
s and p waves (see text). 
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-26.28 


-26.28 


-26.28 
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-79.6 


-33.39 


-79.6 
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3.5 
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33.39 
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u so 


3.5 


3.5 
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3.5 






1.7 


1.7 


2.5 


1.7 



is a non-relativistic potential reproducing proton-proton 
and neutron-proton scattering data for energies from to 
350 MeV, neutron-neutron low-energy scattering data, as 
well as the deuteron properties. 



B. Neutron- 10 Be potential 

For the neutron-core interaction we have constructed 
an ^-dependent potential of the form: 



vW(r) = vM(r) + v£(r)t. 



(5) 



where t is the neutron-core relative orbital angular mo- 
mentum and s n is the spin of the neutron. 

The central (V c (»")) and spin-orbit (Vy (?")) radial 
potentials are assumed to have a gaussian shape. In this 
work four different 10 Be-neutron potentials (labeled as 
/, //, /iT, and IV) will be used. Their parameters for 
£=0, 1, and 2 are given in table HI and they are adjusted 
to reproduce the spectrum of n Be. Contributions from 
partial waves with £ > 2 are negligibly small and not 
included. 

Potential I is built as follows: the range of the inter- 
actions is taken equal to 3.5 fm, that is the sum of the 
rms radius of the core and the radius of the neutron. For 
s-waves the strength is hxed to fit the experimental neu- 
tron separation energy of the l/2 + -state in n Be (—0.504 
MeV [43| )• For p- waves the two free parameters (central 



and spin-orbit strengths) are adjusted to reproduce the 
experimental neutron separation energy of the l/2 _ -state 
in 11 Be (—0.184 MeV j43|), and simultaneously push up 
the 3/2~ state, which is forbidden by the Pauli principle, 
since it is occupied by the four neutrons in the 10 Be core. 

For the d-states it is well established that n Be has a 
5/2 + resonance at 1.28 MeV (energy above threshold) 
[4J, [43, |46[ . The strength of the <i5/2-potential is then 
fixed to reproduce this resonance energy (as a pole of 
the S- matrix), leading to a resonance width of 0.4 MeV. 
The most likely candidate as spin-orbit partner of the 
5/2+ state is the known 3/2 + -resonance at 2.90 MeV 
(above threshold) [40]. A gaussian with a range of 3.5 
fm fitting such 3/2 + energy is giving rise to a very broad 
resonance of roughly 1.5 MeV. Experimentally, the 5/2 + 
and 3/2+ states at 1.28 MeV and 2.90 MeV are rather 
narrow (about 100 keV) (43[. For this reason we have 
reduced the range of the d.3/2 neutron-core interaction to 
1.7 fm, such that the width of the 3/2+ state at 2.90 
MeV is also 0.4 MeV. These conditions lead to a central 
(Vc ~ ) and spin-orbit (VJ ~ ) radial potentials made 
as a sum of two gaussians, whose strengths and ranges 
are given at the bottom of table HI 

In principle the computed widths of the 5/2 + and 
3/2 + resonances can be reduced by simply using smaller 
ranges for the corresponding gaussians. However, to ob- 
tain widths similar to the experimental ones, unrealistic 
ranges are needed, and we have then preferred to use 
a d-wave interaction with ranges similar to the ones for 
the s and p potentials. This disagreement indicates that 
these n Be states, beside the dominating single-neutron 
d-waves, have admixtures of the 10 Be core excited 2 + 
coupled to the single-neutron s-wave. 

Among the different partial wave neutron-core poten- 
tials described above, the s-wave potential is probably 
the most crucial, since it determines the properties of 
the n Be ground state. The s-wave neutron- 10 Be compo- 
nents are expected to give a large contribution not only 
to the 12 Be ground state but also to most of the excited 
states. To test the dependence of the results on the de- 
tails of this potential, we have in potential II increased 
the range of the s-wave neutron-core interaction up to 
4.5 fm, modifying the strength to keep the energy of the 
1/2+ state in n Be at the experimental value. 

When constructing potential I the range of the d 3 / 2 
potential was reduced to obtain a width for the 3/2+ 
resonance in better agreement with the experiment. To 
test the importance of this choice, we have in potential 
III increased the range of the d^/2 interaction to 2.5 fm. 
The width of the 3/2 + resonance (at 2.90 MeV above 
threshold) is then 1.1 MeV. 

The 10 Be core has two neutrons occupying the S1/2- 
shell and four neutrons in the p3/2-shell. Therefore, the 
neutron-core interaction can not bind the neutron into 
one of these states due to the Pauli principle. To treat 
this fact in a more careful way, we have constructed a 
new neutron-core potential such that there is a deeply 
bound Si/2-state (forbidden by the Pauli principle) and 



TABLE II: Strengths (S) and ranges (6) used for the gaussians 
describing the deep neutron- 10 Be interaction holding Pauli 
forbidden states for the S1/2 and p 3 / 2 -waves. The last column 
gives the single nucleon rms radius for each wave. 



S (MeV) 



b (fm) 



(rl) 1/2 (fm) 



s i/2 _wa ve 
p 3/2 -wave 



-139.84 
-75.5 



1.78 
2.3 



1.26 
2.72 



a second bound s^-state at the experimentally known 
energy of —0.504 MeV. In the same way the p 3 / 2 interac- 
tion has been made such that there is a bound state (also 
forbidden by the Pauli principle) at —6.81 MeV, that is 
the experimental neutron separation energy in 10 Be [4 71 ]. 
These two potentials are also taken to be gaussians, and 
their strengths and ranges are given in the second and 
third columns of table |TTJ The last column gives the 
rms radii for a neutron sitting in the Sj / 2 -shell or in the 
P3/2-shell in 10 Be. The parameters used for the gaussian 
potentials have been adjusted such that, together with 
the proper S1/2 and P3/2 energies, the computed single 
nucleon rms radii give rise to charge and mass rms radii 
for 10 Be in agreement with the experimental values, i.e. 
2.24 fm and 2.30 fm, respectively [48J |. 

We have then constructed potential IV using the S1/2 
and p 3 /2 interactions as the phase equivalent potentials 
[49j of the ones described above and given in table [Til 
These potentials have exactly the same phase shifts for all 
the energies as the original ones, but the Pauli forbidden 
states have been removed from the two-body spectrum. 
Potential IV is completed with the same P1/2 and d- 
wave interactions as in potential /. In particular, the 
central and spin-orbit parts for p-waves in potential IV 



r(«=l) 



(2V t 



P3/2 



V Pl/2 )/3 and V s 'o 



(£=1) 



V Pl/2 )/3, where V P3/2 and V Pl/2 are the p 3 / 2 and 



are given by: V c 

2(V P3/2 

P1/2 potentials as described above 

The different 10 Be-neutron interactions described in 
this section reproduce reasonably well the known spec- 
trum of n Be up to an excitation energy of 3.41 MeV. 
However, it is well established that n Be has a 3/2~ res- 
onance at 2.19 MeV (above threshold) [50, [5J, [53] , which 
has been ignored in our analysis. This is because a 10 Be 
core in the + ground state can not produce a low-lying 
3/2 _ state in 11 Be. The first allowed p 3 / 2 -shell where the 
halo neutron could sit is too high (even if a large 10 Be 
deformation was assumed). In fact, the 3/2~ state in 
11 Be very likely corresponds to 9 Be (whose ground state 
is 3/2~) and two neutrons in the sd-shell [44j. In other 
words, one of the neutrons in the fully occupied p 3 / 2 -shell 
in 10 Be has to jump into the sd-shell, which means that 
a description of this 3/2~ resonance as a 10 Be core plus 
a neutron requires the core in a negative parity excited 
state. Another possibility is to have the 10 Be core in 
the excited 2 + state and the remaining neutron in the 
Pi/2-shell. Therefore, when investigating how an inert 
core three-body model describes the 12 Be properties this 



TABLE III: Components included in the calculations for the 
+ states. The upper part corresponds to the first Jacobi set 
(x between the two neutrons). The lower part corresponds 
to the second and third Jacobi sets (x from core to neutron). 
The 6' column gives the maximum value of the hypermo- 
mentum used for each component. The last five columns give 
the contribution from each component to the + wave func- 
tion for potentials I, II, III, IV, and I+i>s, respectively. The 
two numbers for each component correspond to the 0^ and 
0^~ states, respectively. 



TABLE IV: The same as table [TTT] for the 2+ state in 12 Be. 
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90.2 
48.3 

8.2 



88.5 
52.5 
10.3 



2 2 



82 



50.2 46.0 
1.6 1.2 
1.5 1.6 



88.9 
49.3 
9.4 
49.2 
1.8 
1.6 



88.7 
53.9 
9.4 
43.5 
1.9 
2.7 
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1/2 
1/2 
1/2 
1/2 
1/2 



118 



80 



80 



82 



82 



75.5 
14.5 
6.4 
29.8 
7.1 
48.1 
9.7 
4.3 
1.2 
3.3 



71.5 
23.3 

7.9 
26.5 
9.0 
44.3 
10.2 
3.0 
1.4 
2.8 



70.1 
15.6 
7.7 
29.2 
8.7 
46.9 
11.9 
4.8 
0.7 
3.4 



66.6 
22.0 
11.1 
29.8 
7.9 
41.6 
12.6 
3.6 
1.6 
2.9 



t'8 



86.5 
49.1 
11.9 
49.5 
1.6 
1.5 



68.9 
16.3 
9.0 
29.5 
10.9 
48.5 
10.0 
3.6 
1.3 
2.2 



two-body 11 Be state has to be excluded. We shall later 
discuss the consequences of this exclusion. 
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0.2 
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0.1 
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38.6 39.4 
8.6 6.9 
37.1 37.9 



8.5 
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6.7 
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1.3 
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0.2 
5.3 
0.02 
0.1 
44.2 
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42.8 
2.9 
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0.2 
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0.01 
0.1 
43.8 
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41.1 
2.4 
0.3 



VH 



34.2 
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9.1 
55.3 
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1.3 
0.2 
0.3 
3.9 
0.1 
0.3 
42.1 
5.7 
40.1 
5.5 
0.1 



TABLE V: 
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IV. STRUCTURE OF 12 BE 

To solve the angular part of the Faddeev equations 
(Eq.©) it is necessary to specify the components in- 
cluded in the hyperspherical harmonic expansion of the 
angular eigenvalues. These components should be con- 
sistent with the total angular momentum and parity of 
the system. For 12 Be two + , one 2 + , and one l - bound 
states are experimentally known. We find all of them and 
predict the existence of an isomeric 0~ state [53[ . The 
components included for them in the numerical calcula- 
tions are given in the first five columns in tables IIII1 IIV1 
and IVIi respectively. The upper part of the tables refer 
to the components in the first Jacobi set (x between the 
two neutrons) , while the lower part gives the components 
in the second and third Jacobi sets (x from the core to 
one of the neutrons). 

The quantum numbers £ x and £ y are the orbital an- 
gular momenta associated to the x and y Jacobi coordi- 
nates, and they couple to the total orbital angular mo- 
mentum L. The spins of the two particles connected by 
the x coordinate couple to s x , that in turn couples with 



the spin of the third particle to the total spin S. Finally 
L and S couple to the total angular momentum of the 
three-body system. An additional quantum number to 
be considered is the hypermomentum (K — 2n + £ x + £ y ) 
whose maximum value (K max ) for each component is cru- 
cial to guarantee convergence of the eigenvalues (A„(p)) 
up to distances where the radial wave functions are neg- 
ligible. The values of K max used in our calculations are 
given by the sixth column in the tables. 



A. Effective potentials 

The angular eigenvalues obtained from Eq. ^ enter in 
the coupled set of radial equations ([3]) as a part of the ef- 
fective potentials V^lifi) = jp (A„(p) + x)- * n Fig 111 we 
show the three most contributing effective potentials for 
the + (upper- left), 2 + (upper-right), 1~ (lower-left) and 
0~ (lower-right) states in 12 Be. The results using poten- 
tial / for the neutron-core interaction and the gaussian 
neutron-neutron potential are shown by the solid curves. 
When potential IV is used (dashed curves) these poten- 
tials are noticeably different in amplitude from the ones 



TABLE VI: The same as table HIT] for the 1" state in 12 Be. 
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with potential /, but, in general, have similar shape. The 
following effective adiabatic potentials behave very sim- 
ilar in both cases. When the neutron-core potentials II 
and III are used, only minor differences in the adiabatic 
potentials are found compared to the ones obtained with 
potential 7. All the neutron-core interactions give rise to 
indistinguishable effective potentials at large distances. 
For the + and 0~ cases the deepest and second deepest 
effective potentials go asymptotically to —0.504 MeV and 
—0.184 MeV, respectively, which are the two-body bind- 
ing energies of the bound states in 11 Bc. For the 1~ and 
2+ the two deepest effective potentials go asymptotically 
to -0.504 MeV and the third deepest to -0.184 MeV. 



B. Three-body energy spectrum 

After solving the coupled set of radial equations 
(Eq.Q) we obtain a series of 12 Be bound states whose 
two-neutron separation energies are given in table IVIII 
Columns 2 to 5 show the results obtained with the 
10 Be-neutron interactions given in table [IJ The gaussian 
neutron-neutron potential is used. In order to test the 
role played by the short-distance details of the nucleon- 
nucleon interaction, we give in the sixth column the spec- 
trum obtained when potential I is combined with the 
Argonne vg neutron-neutron interaction. The last col- 
umn gives the available experimental values. These re- 
sults correspond to calculations without fine tuning the 
effective three-body potentials with a three-body force 
{V 3 b(p)=0 in Eq.©). 

The bound states found can be understood assuming 
a simple extreme single particle model: The ground + 
state should correspond to a configuration where the two 
neutrons occupy the two S1/2 single-neutron states. The 
1~ and the 0~ can be understood as one neutron in the 
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FIG. 1: Three lowest Adiabatic effective potentials for the 
+ , 2 + , and 1~ states in 12 Be. Results using the neutron- 
core potentials I (solid curves) and IV (dashed curves) are 
shown. 



TABLE VII: Spectrum of 12 Be for the different neutron-core 
interactions (see tableUJ. The numerical results have been ob- 
tained without inclusion of a three-body potential in Eq.©. 
In the sixth column potential I has been used together with 
the ug-Argonne nucleon-nucleon interaction. The energies are 
given in MeV. 
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S]/2 wave and the other one in a P1/2 state. The 2 + 
state appears with one neutron in the S1/2 wave and the 
second in the <7 5 / 2 state, and finally, the second + state 
corresponds to both neutrons in p!/ 2 -waves. 

For the Of ground state all the computed energies are 
similar, except for potential 77, that is clearly less bound 
than the rest. In potential 77 the s-wave neutron-core 
potential employed a larger range (see table H]). Since the 
structure of the 0^~ state is expected to be dominated 
by the s-wave components, an increase of the range for 
this wave increases the spatial extension of the system 
and subsequently decreases the binding energy. In po- 
tential IV the s-wave interaction has been constructed 
in a completely independent way compared to potentials 
7, 77, and III. The fact that in this case the binding en- 
ergy of the ground state is similar (actually slightly more 



bound) than for potentials / and III, indicates that a 
range of 3.5 fm for the s-wave interaction could be more 
appropriate. 

The p-wave potential is identical for the cases /, II, 
and III. Potential IV also uses the same P\j 2 interac- 
tion, but different p 3 / 2 , which however is expected to play 
a minor role due to the Pauli principle. This explains 
why the energy of the 2 state is very much the same 
for all neutron-core potentials, since this state should be 
dominated by a pp configuration. 

For the 1~ (or 0~) and 2 + states a neutron in the 
S]/2 wave is combined with the second neutron either in 
the P1/2 or the d 5 / 2 state, respectively. These states are 
less bound than the ground state. The effect of this is 
more important than the confining effects of the p and 
d centrifugal barriers, and consequently the 1" (or 0~) 
and 2 + states are more extended than the ground state. 
Therefore the effect of the larger range in potential II 
for the s-wave interaction should be smaller. Actually, 
the 1~ energy is the same for all the three potentials. 
Only potential III, for which the d^/ 2 wave was modified, 
produces a slightly more bound 2 + state. 

In potential IV the d-wave interaction is the same as 
in potential /. Also the s-wave and p-wave potentials, 
although obtained in a completely different way, produce 
similar + states. However, the 1~ state, and especially 
the 2 + state, are more bound than with potential /. This 
indicates that the sp and sd interferences for these states 
differ between potentials TV and /. When the gaus- 
sian neutron-neutron interaction is substituted by the 
i>8-Argonne potential (sixth column in table IVIip simi- 
lar results are found. Only a little more binding is found 
for the 1~, 0~ and 0^ states. 

In summary, the energies obtained with the differ- 
ent potentials are quite stable, and the small differences 
found for some of the cases are insignificant, especially 
when taking into account the level of accuracy of the cal- 
culations at this stage. In fact, when compared to the 
experimental energies (last column of table IVII[) we see 
that in some cases the computed energy really disagrees 
with the experimental value, and even the ordering of 
the levels does not agree with the experiment. However, 
the ordering of the computed levels agrees with the ex- 
treme single particle model. The ground state should be 
the one having the two neutrons in the deepest single 
neutron level (the Of state having both neutrons in the 
Si/2-wave). The first excited state should appear when 
one neutron jumps into the next single neutron level (the 
1~ or the 0~ state having one neutron in the Si/2-wave 
and the other one in the pi/2-wave). Finally, the states 
corresponding to the p 1 / 2 — P\/ 2 and the S\ m — d 5 / 2 con- 
figurations (Oj and 2 + states) should in principle be less 
bound than the 1~ and the 0~ states. 



C. Angular momentum decomposition 

It is important to keep in mind that part of the con- 
tributions arising from core deformation have been effec- 
tively taken into account by fitting the parameters in the 
two-body potentials to the 11 Be data. In this sense, the 
different states found in the 12 Be spectrum are sensitive 
to such deformation. However, with such two-body in- 
teractions and the core treated as an inert particle, the 
additional contributions arising from the 2 + excited state 
in 10 Be are expected to play a minor role for both the Of 
and 1~ states. For the ground state because two neu- 
trons in the s^-shell and a core in a 2 + state can not 
produce total angular momentum zero. For the 1 _ state 
because, although a core in a 2 + state and the two neu- 
trons in the S1/2 — P1/2 configuration can produce a total 
angular momentum 1, this structure is obviously less fa- 
vorable energetically than with the core in the ground 
state. However, for the 2 + state it could be energetically 
efficient to excite the core into the 2 + state and get some 
extra binding by placing the two neutrons in a S\i 2 — s 1 i 2 
configuration. The same could happen for the second + 
state, that could find it favorable to excite the core into 
the 2 + state and place the neutrons in the S1/2 — d§/ 2 
configuration. Therefore, the disagreement between the 
computed and experimental binding energies of the 2 + 
and 2 states can be a signal of the importance of core 
excitations for these two particular states. 

In any case, the deviations between computed and ex- 
perimental energies are quite common when three-body 
calculations are performed with bare two-body interac- 
tions. All those effects going beyond the two-body corre- 
lations have obviously not been considered in the calcula- 
tions. Among them are those arising from contributions 
of core deformation and/or core excitation. The usual 
cure for this problem is including an effective three-body 
potential (V^bip)) that simulates the neglected effects. 
Since these effects should appear when all the three par- 
ticles are close to each other, the three-body potential 
should be of short-range character. 

In our calculations we have used a gaussian three-body 
force, whose range has been taken equal to 4.25 fm, that 
is the hyperradius corresponding to a 10 Be-core and two 
neutrons touching each other. The strength of the gaus- 
sian is adjusted to match the experimental energies given 
in the last column of table IVII1 

Since the CP-state is unknown experimentally, we can 
not use its energy to adjust the strength of the three 
body force. However, in this case we expect the three- 
body correction to be unimportant due to the similarity 
of the 0~ state to the composition of the l~-state where 
the experimental energy is reproduced without any three- 
body force. This indicates that the computed 0~-energy 
also is close to the correct value. 

Including the three-body potentials we have the four 
potentials specified in table [H and in addition potential 
/ combined with the ws-Argonne neutron-neutron poten- 
tial. We then arrive at the contributions of the different 



TABLE VIII: Root mean square radii (in fin) for the different 
bound states in 12 Be with the four neutron-core interactions. 
A gaussian three-body forced is included to fit the experi- 
mental binding energies. The range of the three-body force 
is 4.25 fm. Sixth column is as potential I + Vs Argonnc 
nucleon-nucleon interaction. 
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partial wave components given in the last five columns 
of tables EI] to EH for the 0+, 2+, 0" and 1~ states. The 
dominating components are in all the cases as expected 
from the extreme single particle picture. When the three- 
body wave functions are written in the second or third Ja- 
cobi set (x from core to neutron, lower part of the tables) 
the dominating components are the £ x =£ y =0 (~ 70%) 
for the Of ground state, the £ x =£ y =l (~ 75%) for the 



Of state, the 



=0, 



=2} and 



=0} components 



(~ 90% in total) for the 2+ state, and the {4=0,4=1} 
and {£ x —l,£ y —0} components (~ 98% and ~ 99% in to- 
tal) for the l - and 0~ states. Although the precise num- 
bers can change a little from one potential to another, 
the differences are not significant. The residual contri- 
butions are more relevant for the + states. The ground 
state has a p-wave contribution of about 13% to 19%, 
and a d-wave contribution of about 10% to 13%. The 
second + state has an s-wave contribution from 15% to 
23%, and a d-w&ve contribution from 6% to 8%. 

In [12j the configuration with the two outer neutrons 
in the p-shell for the + ground state amounts to 32% of 
the wave function, which is clearly larger than the 13%- 
19% obtained in this work. In [2] this value is given to 
range between 20% and 40%. Correspondingly, the con- 
tribution from configurations with the two neutrons in 
the sd-sflell is larger in our calculation than in [2, LL2J- 
In these two references the individual contributions from 
s and d- waves are not given. In general, for both n Li 
and 12 Be the s-wave components are larger in cluster 
models than in shell model calculations. One reason for 
12 Be is that the d-waves are underestimated in the clus- 
ter model due to the neglect of the core-excited 2 + state. 
Another reason could be that large spatial extension is 
harder to describe in shell models than in cluster mod- 
els. Since s-waves for a given energy extend to larger 
distances than d-waves the shell model tend to underes- 
timate the s-wave components. As we shall show later 
the contributions obtained in the present work are con- 
sistent with the measured invariant mass spectrum given 
in El. 



TABLE IX: For the different computed states in 12 Be, and 
the different neutron-core potentials, root mean square dis- 
tances (in fm) (r nn } and (r C n) , where n and c denote 
an external neutron and the core, respectively. A gaussian 
three-body force is included to fit the experimental binding 
energies. 
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FIG. 2: Contour diagram for the probability distribution of 
Of state in 12 Be. The square of the three-body wave function 
is integrated over the directions of the two Jacobi coordinates. 



D. Wave functions 

The rms radii for the computed bound states are given 
in table IVIIII The results are essentially independent 
of the core-neutron potential used. The only available 
experimental value is the one corresponding to the Of 
ground state [48j, for which a good agreement between 
theory and experiment is found. 

The geometry of the different states is reflected in the 
results given in table iLXl where we give the rms distances 
between the two neutrons ((rj^) 1 ' 2 ), and between the 
core and one of the neutrons ((^n) 1 / 2 ) for the different 
cases. The ground state corresponds mainly to the three 
particles placed in the vertices of an equilateral triangle. 
For the excited states, the smaller binding obviously im- 





FIG. 6: The same as Fig[2]for the (T state in 12 Be. 



FIG. 3: The same as Fig0fc>r the 0^" state in 12 Be. 
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FIG. 4: The same as Fig[2]for the 2+ state in 12 Be. 
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FIG. 5: The same as Fig[2]for the 1" state in 12 Be. 



plies a larger distance between particles. But, as seen in 
tablc lIXl the distance between the neutrons grows clearly 
faster than the one between the neutron and the core. In 
fact, for the less bound excited states (the l - and 0~ 
states) the distance between the two neutrons is roughly 
a factor of two larger than the corresponding distance 
in the ground state, while the neutron-core distance in- 
creases by a factor of 1.5. 

The spatial distribution of the three constituents can 
be better seen in Figures O to [SI where we show the prob- 
ability distribution for the 0^~, 0^, 2 + , 1 _ and - states 
in 12 Be. The probability distribution is defined as the 
square of the three-body wave functions multiplied by the 
phase space factors and integrated over the directions of 
the two Jacobi coordinates. In particular, in the figures 
we have chosen the first Jacobi set, such that the distri- 
butions are plotted as function of the distance between 
the two neutrons (r nn ) and the distance between the core 
and the center of mass of the two neutrons (r Cj „ n ), re- 
spectively The 0+ and the 2 + states have the largest 
probabilities when the two neutrons are separated from 
each other almost the same distance as from the core 
(about 2.2 fm and 2.6 fm, respectively). On the other 
hand, 0^", 0~ and 1~ states have the largest probabilities 
when the two neutron are well separated from each other 
(about 5.5 fm) while they are closer to the core (about 
3.5 fm). The 2 + state has also a secondary probability 
maximum with this last geometry. 



V. TRANSITION STRENGTHS 

The bound states can only decay electromagnctically 
The corresponding observable transition probabilities are 
critically depending on the structures. Thus they provide 
experimental tests and we therefore compute the lifetimes 
for future comparison. The selection rules determine the 
dominating transitions which can be of both electric and 
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magnetic origin. 

A. Electric transitions 

The electric multiple operators are defined as: 

A 

M»(EX) = eJ2 z i r i Y ^An) 



(6) 



where A is the number of constituents in the system, each 
of them with charge eZ t , and where ?% is the coordinate 
of each of them relative to the A-body center of mass. 
The electric multipole strength functions are defined 



B(EX,h^I f ) = YsKIfMflM^EXpiMJl 2 

= ^jKIfWMiEXWi)] 2 (7) 

Also, following [401, we consider the monopole transi- 
tion operator: 



M o (E0) = e^TZi 



(8) 



In our three-body model, where A — 2 constituents are 
assumed to form a well defined core, the operators Q 
and © contain one term directly associated to the three- 
body system having the form: 

3 

M M {EX) = eJ2ZirtY x>IM (h) (9) 

i=l 
3 

Mo(EO) = e^Z 4 rf (10) 



where i labels the three constituents. In addition to 
Eqs.© and (jTUJ) there is a contribution from the intrinsic 
core multipole operators. In particular, for the electric 
dipole and quadrupole operators the precise expressions 
are given by Eqs. (|A3|) and (|A4[) . For the quadrupole 
case, Eq. (|A4[) . there is an additional contribution arising 
from the coupling between the electric dipole operator of 
the core and its motion around the three-body center-of- 
mass. 

When an inert core has spin zero, as in the present 12 Be 
model, the contributions from the intrinsic core multi- 
pole operators M^EX, core) are zero. These terms con- 
tribute only when core excitations are included. For 12 Be 
the first excited state of the 10 Be-core is a 2 + state at 
about 3.4 MeV. Its effect on the B{E1) strength must be 
small, since the M. ll {E\,core) operator does not couple 
the + and 2 + core states. Therefore the B(E1) strength 
is expected to be less dependent on the contribution of 
the 2 + excited state in the core than the B(E2) strength, 



where the A / f M (£'2, core) operator does couple the + and 
2 + states in the core. 

Also, the M o(750, core) operator is not coupling the 
+ and 2 + states of the core. However, the expectation 
value of this operator between the + ground state of the 
core is not zero (it is actually related to the rms radius of 
the core). Nevertheless, for an inert core with spin zero, 
the total nuclear wave function factorizes into the three- 
body cluster wave function and the core wave function. 
In this way, the expectation value of Mo(EO,core) be- 
tween wave functions corresponding to different states is 
automatically zero due to the orthogonality of the three- 
body cluster wave functions. 

For n Be ( 10 Be+n) the experimental value of the 
B{El,\~ -> i + ) transition strength is 0.115 ± 0.010 
e 2 fm 2 [54j . Direct application of Eq. ([9]) assuming zero 
charge for the neutrons fails in reproducing this value. 
However, it is well known that to take into account 
the distortion or polarization of the core one has to 
include an effective charge for the neutrons [55|, l56j |- 
Due to the square in Eq.Q two different neutron effec- 
tive charges are found to fit the experimental value for 
B(E1, \~ -> \ + )- Z n = 0.28 and Z n = 0.52. These two 
charges give rise to slightly different mean square charge 
radii for the ground state in n Be of about 2.7 fm and 
3.0 fm, respectively. The 10 Be-neutron potential only 
has marginal effect and only 2.7 fm is_consistent with 
the experimental value of 2.63±0.05 fm [48|. Therefore a 
neutron effective charge of Z„=0.28 has been used in the 
following calculations. 

In table|X]we gi ye the computed B(EX) values for tran- 
sitions between the states in table IVIII and the monopole 
transition matrix element M(E0) = \(0^\M o (E0)\0f )|. 
An effective three-body force has been included to fit 
the experimental binding energies for the different 12 Be 
states. The results are quite stable for the different 
neutron-core potentials used. The only exception are the 
transitions in which the 0J" state is involved. This state is 
the one showing the most important dependence on the 
potential used. In particular potentials 77 and IV pro- 
duce an s-wave content (lower part of table IIII[) about 
50% larger than potentials I and 777. As seen in the 
table the results are sensible to this difference, since the 
values found for £(752,2+ -> 0+), B(E1,0+ -> 1~) and 
A4(E0) with potentials II and IV clearly differ from the 
ones with I and 777. 

For the 73(751) transitions, the experimental value of 
0.051 ± 0.013 e 2 fm 2 corresponding to the 0^ — > 1 _ 
transition [32], agrees well with the results given in the 
first row of table |Xj The agreement is equally good for 
all the neutron-core potentials used. Only for potential 
IV the computed value is a bit higher compared to the 
other cases, but the computed result is still lying within 
the experimental error. Therefore the three-body model 
used reproduces well the experimental value, although it 
is worth emphasizing that the effective neutron charge 
also is introduced to simulate neglected effects of core 
deformation and polarization. The small effect from the 
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TABLE X: B{E1) (in e 2 fm 2 ), B{E2) (in e 2 fm 4 ) transition strengths between the computed states in 12 Be given in table I VIII 
and M(E0) (in e f m ) between the two + states. An effective three-body force has been used to fit the experimental binding 
energies. 
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0.016 


0.010 


0.015 


— 


B(E1,2+ 


-n 


0.0054 


0.0074 


0.0079 


0.018 


0.0081 


— 


B(E2,2 + 


-0+) 


3.14 


3.50 


3.52 


3.93 


3.04 


5.5-8.2 (i,) 


B{E2,2 + 


-0+) 


0.024 


0.32 


0.044 


0.12 


0.047 


1.40±0.12 (c) 


M(E0,0+ 


-o+) 


1.89 


0.60 


1.77 


1.13 


1.99 


0.87±0.03 (c) 



'"' from [3^| . ' estimation from experimental data in [3j|, ' from [4i 



three-body potential can be attributed to the choice of 
minimal structure, which means no dependence on an- 
gular momentum quantum numbers. The only essential 
effect is an adjustment of the three-body energy which 
in turn may have an effect on the spatial extension, but 
leaving all structures unchanged. For the other two El 
transitions experimental data are not available. 

Using the same effective neutron charge we get the 
B(E2) values in the central part of table |X] Assuming 
that the 2 + state in 12 Be is a rotational state built on the 
ground state, an experimental value of 2.00±0.23 is given 
in [31| for the deformation length 5 in 12 Bc. This assump- 
tion means that the B(E2, 2 + — > Of ) value is given by 

\ (^ZeS) 2 R 2 , which equals 5.5 ± 1.3 e 2 fm 4 or 8.2 ±1.9 
e 2 fm 4 for R « 1.2A 1 / 3 fm or R = ^^(r 2 ) 1 / 2 , re- 
spectively. Here (r 2 ) 1 ' 2 is the experimental charge r.m.s. 
radius of 12 Be. Both values are higher than obtained 
numerically (first line in the central part of the table) . 

Also the recently measured value of 1.40 ± 0.12 e 2 fm 4 
[ioj for B(E2, 2+ -> 0j) is much larger than all the com- 
puted results in the table, which varies about one order 
of magnitude with the different neutron-core potentials. 
Furthermore we see that the computed B(E2, 2+^0^) 
is much larger than B(E2,2 + — * 0~£). This arises from 
the fact that the 2 + state is dominated by the sd interfer- 
ences in the second and third Jacobi sets (lower part of ta- 
ble HY|, while the O2" state is dominated by a pp configu- 



ration (lower part of tablc lIlT]) . When the core is infinitely 
heavy these dominating configurations must give a van- 
ishing contribution from Eq.@ to the B(E2,2 + -> Of) 
value. Thus, the values in table |X] have to be small and 
very sensitive to the non-dominating components in the 
two states. 

Small variations in the contribution of some of these 
smaller components can produce large relative changes in 
the computed transition strength. In particular, for the 
Oj state, the £ x — 1 components clearly dominate but 
a substantial probability appears in the {£ x — 0, S = 0} 
component, which coincides with a large probability for 
the {£ x = 0, S = 0} component in the 2 + state, see the 
lower part of table IIVI For potentials II and IV this s- 
wave component in the Of state (table HTl)) is significantly 



larger than for the other three potentials. This implies 
a larger overlap, and therefore a larger B(E2) transition 
strength. 



At this point it is important to emphasize that the 
computed results are obtained for an inert core. For 
i?2-transitions a non-negligible contribution from the 2 + 
excited state in the 10 Be-core is expected. The second 
term in the electric quadrupole operator (|A4[) , together 
with the existence of the core excited 2 + -state, give rise 
to another non-vanishing contribution. In fact, for 10 Be 
the experimental B(E2,2 + — > + ) transition strength is 
known to be 10.4 e 2 fm 4 [53]. Thus the computed three- 
body B(E2) values should be supplemented by the con- 
tribution of 10.4 e 2 fm 4 multiplied by the weight factor 
corresponding to the admixtures of the 2 + core excita- 
tion in the 12 Be states. In fact, this contribution should 
provide most of the strength for the 2 + — ► 0J" transi- 
tion. In contrast, the contribution from the last term in 
the quadrupole transition operator (|A4[) is expected to 
be small, since M f _ l (El,core) cannot couple + and 2 + 
states. 



In [40j an experimental value of 0.87 ± 0.03 e fm 2 is 
given for the monopole transition matrix element be- 
tween the two + states. This observable gives informa- 
tion about the relative structures of the these two states. 
The computed results are given in the last row of tablelXl 
and they are clearly higher (except for potential II) than 
the experimental value, but similar to the value of 1.7 e 
fm 2 obtained in Ref. [5g|. However, as seen in the table, 
the computed results are very sensitive to the contribu- 
tion of the relative neutron- 10 Be s-wave. An increase 
from about 15% to 23% reduces the computed A4(E0) 
by a factor of 3 (going in fact through the experimental 
value). This s-wave contribution could easily change sig- 
nificantly when core excitations are included in the Of 
wave function, for which, as mentioned at the beginning 
of subsection II V Cl such excitations could play a relevant 
role. 
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B. Magnetic transitions 

For the unobserved predicted 0~ state the dominat- 
ing multipolarities for its possible decays are of magnetic 
character. Only Ml and Ml transitions to the 1~ or the 
2+ states are possible. 

The magnetic multiple operator is defined as: 



M^MX) 



eh 
2Mc 

A + l 



v/A(2A + l)5y, 



,A-1 



(Yx-is) + 



2.9, 



(i) 



A + l 



(Xx-ij) 



(11) 

(A-l,l)Aju 



where the constants g s and ge depend on the constituent 
particles i. The magnetic multipole strength functions 
are defined as for the electric case (see Eq. ((?])). 

In particular, the .Ml and A42 operators involved in 
the magnetic dipole and quadrupole decay of the 0~ state 
into the 1~ and 2 + states are 



M^Ml) = 



M»{M2) = — -= 




gfl). 



(12) 



Yi tV (SU) 



(13) 



where q labels the spherical component of an operator. 

We can identify a source of uncertainties in the tran- 
sition strength estimates, which comes from the effective 
values of the g- factor in these expressions. 

The core has angular momentum zero and thcrc- 

{ c) (c) 

fore a vanishing effective spin g s -factor and g) = 4. 
We also use the free value of g s =— 3.82 and we use 
again an effective neutron charge g £ n — 0.28. The tran- 
sition operators arc then defined and we can compute 
the B(M1,0~ -> l - ) and B(M2,0 _ -> 2+) transition 
strengths. The results are given in table IXI[ where we 
observe that the computed values are rather independent 
of the particular neutron-core potential used. Only for 
potential IV the results deviate by up to factor of two 
from the other estimates. This is reflecting the fact that 
potential IV is producing a different distribution of the 
weights between the components in the l - and 2 + states 
(see tables IVHarid lIV|) . 

As mentioned above, the effective values of the g- 
factors are rather uncertain and spin polarization could 
reduce g s by a factor of 2, change g[' by perhaps 10 %, 

and vary g^ n from the assumed effective value (see also 
the discussion of the empirical evidence in [5!|). As dis- 
cussed in [531 ] , the computed magnetic strength has a very 
limited dependence on the precise values used for the or- 
bital gyromagnetic factors ge, and they are mainly only 
sensitive to g s . In particular, use of gs =— 2.0 reduces 
the computed magnetic strengths by a factor of about 3. 



TABLE XI: B(M1) (in e 2 fm 2 ) and B{M2) (in e 2 fm 4 ) transi- 
tion strengths between the computed states in 12 Be given in 
table rvTTl An effective three-body force has been used to fit 
the experimental binding energies. 



Wi 



Wh Win Wiv Wi+ 



VH 



B(M1,0- 

B(M2,0~ 



1-) 

2+) 



0.060 0.061 
0.94 0.78 



0.059 
0.78 



0.049 0.061 
0.43 0.85 



In summary, the decay possibilities for the 0~ state are 
limited to magnetic transitions. Even with fairly reliable 
estimates for the transition probabilities it is not easy to 
give an accurate prediction for the lifetime of this state. 
If 0~ is below the 1~ state only decay into 2 + is possi- 
ble and the resulting lifetime can be estimated to be of 
the order 10 -8 sec. However, if 0~ can decay into the 
1~ state the possibly very small energy difference causes 
a large uncertainty which could reduce the lifetime to 
about 10 -11 sec or perhaps even somewhat smaller. In 
any case these estimates justify the classification of this 
new 0~ state as an isomer in 12 Be. 



VI. INVARIANT MASS SPECTRUM 

One of the tools for studying the structure of three- 
body halo nuclei was breakup reactions on a target. By 
relatively fast removal of one of the constituents the re- 
maining two are essentially left undisturbed and mea- 
surements of their momenta then provide fairly direct 
information about the initial state. We use this tech- 
nique to investigate the two-body substructures of 12 Be. 
The experimental relative 10 Be-neutron energy spectrum 
(or invariant mass spectrum) after 12 Be breakup on a 
carbon target at a beam energy of 39.3 MeV/nucleon is 
shown in ref.[13J. Both fragments, 10 Be and neutron, are 
simultaneously detected, and their relative decay energy 
is reconstructed from the measured momenta. The final 
10 Be-neutron states are then restricted to unbound states 
in n Be. 

This spectrum can easily be computed by use of the 
sudden approximation. We assume that the target trans- 
fers a (relatively large) momentum to one of the parti- 
cles (one of the neutrons) in the projectile which then is 
instantaneously removed. The remaining two particles, 
the 10 Be-core and the second neutron, are just specta- 
tors in the reaction, meaning that they continue their 
motion undisturbed without further interaction with the 
removed particle. They do, however, continue to interact 
with each other. 

Under this assumption the differential cross section of 
the process takes the form: 



d 6 a 



E 

M < 



E 



(e ik «- y X° s 



'U' 



■(k x ,x)\V JM ( X ,y)) 



(14) 



where ty JM is the three-body wave function with total 
angular momentum J and projection M , k x and k y are 
the momenta associated to the Jacobi coordinates x and 
y, s x and s y are the coupled spin of the spectator two- 
body system and the spin of the removed particle, re- 
spectively, and a x and a y are the corresponding spin pro- 
jections. Finally, w"* is the continuum two-body wave 
function of the spectators ( 10 Be and neutron) in the final 
state. This two-body wave function is computed with 
the corresponding two-body interaction and the bound- 
ary condition at small distance determined to be pre- 
cisely the 10 Be-neutron structure left by removing the 
other neutron in 12 Be, see [601 f° r a detailed description. 
After integration of Eq. (fl"4j) over k y and the angles 
describing the direction of k x we obtain the differential 
cross section da/dk x , which is related to the invariant 
mass spectrum as shown in [36[: 



2000 



da _ E c E n m(M c + M n ) 1 da 
dE nc ~ E c + E n M c M n V x ~dk x 



(15) 



where E c and M c are the core energy and mass, E n 
and M n the neutron energy and mass, E nc the relative 
neutron-core energy, and m is the normalization mass 
used to define the Jacobi coordinates. 

The presence in Eq. (JT5J) of the two-body wave function 
u>J x (k x ,x) makes it evident that the invariant mass spec- 
trum (|15p crucially depends on the final state two-body 
interaction. It is important to note that Eq. (fT4"|) is not 
assuming that the two spectators populate two-body res- 
onances in the final state. The two-body wave function 
u>3 X contributes for any value of the relative two-body 
momentum k x , not only for those where k x matches a 
two-body resonance energy. Except for those very precise 
values of k x , w^* is just an ordinary continuum two-body 
wave function. 

In Fig|7]we show the relative core-neutron energy spec- 
trum after removal of one neutron from 12 Bc on a light 
target according to Eq. ([13)1 . The results for the different 
10 Be-neutron interactions used in this work are shown. 
The computed curves have been scaled to the experimen- 
tal data [13 |. Also, the experimental energy resolution 
has been taken into account by convoluting the computed 
distributions with a gaussian having a full width at half 
maximum (FWHM) equal to OAE 1 / 2 (E in MeV) [jj]. 

As a general result, all the potentials reproduce equally 
well the peak at about 1.2 MeV, and the tail of the dis- 
tribution. In the same way all of them underestimate 
the spectrum in the region around 2.2 MeV and do not 
reproduce the very narrow peak at very low energies. Po- 
tential IV (dot-dashed curve) is reproducing the peak at 
1.2 MeV best of all, but on the other hand gives a worse 
agreement with the experiment in the region between 3 
and 4 MeV. Compared to the results with potential I 
(solid curve), inclusion of the Argonne nucleon-nucleon 
potential (dotted curve) slightly improves the behaviour 
at small energies but also slightly spoils the agreement 
with the experiment at large energies. 



• Exp. [12] 

Potential I 

Potential II 

Potential III 

Potential IV 

Potential I + V t 




3 4 

E . (MeV) 



FIG. 7: Relative energy spectrum of 10 Be+neutron after frag- 
mentation of 12 Be on a light target. The solid, short-dashed, 
long-dashed, dot-dashed, and dotted curves correspond to cal- 
culations using potentials I, II, III, IV, and I + Vs, respec- 
tively. The experimental data correspond to fragmentation of 
12 Be on carbon at a beam energy of 39.3 MeV/nucleon [la ]. 



In Fig|5]we show, for different 10 Be-neutron potentials, 
the contribution from s-, p-, and d- waves to the total in- 
variant mass spectrum. Here s, p, and d refer to the value 
of the 10 Be-neutron relative orbital angular momentum. 
In the figure the dotted and dot-dashed curves give the 
contribution from the d 3 / 2 and d 5 / 2 -waves, respectively. 
We notice that the main peak in the distribution is pro- 
duced by the ds/2-wave. This can be traced back to the 
5/2+ resonance in 11 Be at 1.28 MeV above threshold. 
The 3/2 + resonance, with an energy of 2.90 MeV, helps 
to match the experimental tail of the distribution. The 
s- and p-waves contribute mainly at small energies, and 
they are almost entirely responsible for the distribution 
below 1 MeV. 

From FigJ5]it is now easy to understand that the main 
disagreement between the computed curves and the ex- 
perimental data is due to the absence in the calculation of 
the 3/2 _ resonance in 11 Be, whose energy above thresh- 
old is 2.19 MeV, precisely the energy region where the 
disagreement is found. However, as already mentioned, 
inclusion of this resonance necessarily requires informa- 
tion beyond the present 12 Be three-body model with an 
inert 10 Be core. If we still insist on the picture of a three- 
body system with a 10 Be core, then a 3/2~ state in n Be 
needs either a neutron in the p-shell coupled to the 2 + 
state of 10 Be, or a neutron in the sd-shell coupled to 
10 Be in a negative parity state. Another possibility is of 
course to use a different cluster description for 12 Be, like 
for instance 9 Be in the 3/2~ ground state plus three neu- 
trons. Therefore, except for the lack of this state in the 
11 Be spectrum, the contributions obtained in this work 
for the + ground state are consistent with the experi- 
mental invariant mass spectrum given in [131 ]. 

The highest n Be resonance that has been considered in 
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FIG. 8: Relative energy spectrum of Be+neutron after frag- 
mentation of Be on a light target (thick solid curve), and 
the contributions to it from s- waves (long-dashed), p- waves 
(short-dashed), and d- waves (thin solid). The d-wave con- 
tribution is split into the 1I3/2 contribution (dotted) and the 
d 5 /2 contribution (dot-dashed). The different panels show the 
spectra corresponding to the potentials in tableland the case 
with potential I plus the Argonne v% nucleon-nucleon poten- 
tial. Experimental data from [13 |. 



our calculations is the 3/2 + state at 2.90 MeV. The next 
excited states in 11 Be are the 5/2~ and 3/2~ resonances 
with energies above threshold at 3.49 MeV and 3.46 MeV, 
respectively 52J. For the same reason as the 3/2~ state 
at 2.19 MeV, these two states can not be included in 
our three-body picture with an inert 10 Be core. In any 
case, decay of these two resonances into the ground state 
of 10 Be plus a neutron would give some small contribu- 
tion to the tail of the invariant mass spectrum. However, 
these two states can also decay into 10 Be plus a neutron 
but leaving 10 Be in its 2 + excited state at 3.37 MeV. 
Therefore very little energy is left for the neutron after 
such decay (less than 100 keV 52]), giving rise to the 
experimental sharp peak at very low energies. Theoret- 
ical calculation of this peak requires then also inclusion 
of core excitations in the model. 

In [HI the experimental FWHM of the 10 Be longi- 
tudinal momentum distribution after fragmentation of 
12 Be on a carbon target (beam energy equal to 56.8 
MeV/nucleon) is found to be 194 ± 9 MeV/c. Within 



the sudden approximation, such distribution can also 
be computed from Eq. (fT4")) . see [60] for details. The 
widths of the longitudinal core momentum distributions 
are computed to be 162 MeV/c, 159 MeV/c, 169 MeV/c, 
171 MeV/c, and 166 MeV/c for the neutron-core poten- 
tials J, II, III, IV, and I+Vg, respectively. 

Inclusion of core excitations should help to reduce the 
discrepancy between the computed widths and the ex- 
perimental value. The components having spin of the 
core equal to 2, and coupling to the total angular mo- 
mentum zero of the 12 Be ground state, must necessarily 
contain non-zero values for l x or/and £ y . This means ad- 
ditional centrifugal barriers which attempt to confine the 
three-body system to smaller spatial extension, and sub- 
sequently producing broader momentum distributions. 



VII. SUMMARY AND CONCLUSIONS 

The properties of 12 Be have been investigated assum- 
ing a three-body structure with an inert 10 Be core and 
two neutrons. A description of the nucleus by use of 
the hyperspherical adiabatic expansion method is partic- 
ularly appropriate in this case, since two of the two-body 
subsystems have bound states. The symmetric treatment 
of all the two-body interactions makes it easier to repro- 
duce the correct two-body asymptotics. This could be 
especially important for the weakly bound excited states. 

We constructed four different neutron- 10 Be interac- 
tions, each of them reproducing the known spectrum of 
n Be up to an excitation energy of 3.41MeV. However, 
the 3/2~ resonance in 11 Be has been excluded, since it 
requires the 10 Be core to be in an excited state. This case 
goes beyond the three-body model with an inert core that 
is used in this work. Nevertheless, part of the effects aris- 
ing from core deformation are effectively taken into ac- 
count through the fitting procedure of the neutron- 10 Be 
potential. For the neutron-neutron interaction two dif- 
ferent potentials have been used, both reproducing low- 
energy scattering data. 

We have found two + , one 2 + , one 1~ and one 0~ 
bound states. The first four are known experimentally, 
but not the _ state. The computed binding energies 
agree reasonably well with the experimental value for the 
ground + state and the excited 1~ state. For the second 
+ and the 2 + states a larger discrepancy is found, which 
can be attributed to contributions from core excitations 
additional to the ones masked in the neutron- 10 Be poten- 
tial. To reproduce the experimental two-neutron separa- 
tion energies for these states we included a three-body 
force. 

The dominating components of the wave functions cor- 
respond to the ones expected from the extreme single par- 
ticle model. The ratio of the computed neutron-neutron 
and core-neutron root mean square distances increases 
with the excitation energy. All possible electric and mag- 
netic transition strengths have also been computed. An 
effective charge for the neutron is used to take into ac- 
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count distortion or polarization of the core. This effective 
charge has been obtained by adjusting the calculation to 
reproduce the experimental B(E\) transition strength in 
11 Be, as well as its charge root mean square radius. 

Core excitations can contribute to transition strengths 
in two different ways. There could be a direct contribu- 
tion arriving from the transition between two core states 
or a contribution to the wave function of a 12 Be state 
then leading to an indirect effect on transition strengths. 

The direct contribution of core excitations in the 12 Be 
monopole and dipole strength functions should be negli- 
gible, since r 2 , M.^El) and M^Ml) operators cannot 
couple the J—0 and J=2 states. On the other hand, 
quadrupole transition operators can couple these states 
and there could be an important effect from the excita- 
tion of the core. The possible contribution of core excita- 
tions to wave functions of 2 + and 0^~ states could lead to 
some uncertainties in the calculated transition strengths 
involving those states. Our calculations reproduce rather 
well the available experimental data. 

The relative 10 Be-neutron energy spectrum after 12 Be 
breakup on a carbon target at a beam energy of 39.3 
MeV/nucleon has been computed within the sudden ap- 
proximation. For all the potentials the experimental 
spectrum is rather well reproduced. Only in the region 
around 2.2 MeV the computed curves underestimate the 
experiment. This is due to the absence of the 3/2~ reso- 
nance in our n Be spectrum, which originates from 10 Be- 
core excited states. 

In summary, a frozen-core three-body model is able to 
reproduce most of the properties of 12 Be: ground and ex- 
cited bound states, electromagnetic transition strengths, 
and invariant mass spectrum after high-energy breakup. 
Core excitations, however, are needed to improve the 
two-neutron separation energies for the 0^" and 2 + states, 
to get a better estimate of the quadrupole transition 
strengths and monopole transition matrix element, and 
to fine tune the agreement with the experimental invari- 
ant mass spectrum. 
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FIG. 9: Coordinates used to calculate the Mfj,(E\) operator 
in 12 Be three-body cluster, n, r2 and r c give the positions 
of the three components relative to the center-of-mass of the 
system. The coordinate of the particle i inside the core is 
written as n = r c + r[ where r' t is its position relative to the 
core center-of-mass. 



where Zi is the charge (in units of e) of constituent i, and 
Ti is its position from the A-body center-of-mass. 

For the particular case of 12 Be, we are assuming that 
the system is clustered, showing a three-body structure 
made by a core (containing A — 2 nucleons) plus two 
additional nucleons outside the core. It is then convenient 
to write the coordinates of the particles inside the core 
as: Ti = r c + r^, where r c gives the position of the core 
center of mass relative to the A-body center-of-mass, and 
r'i is the coordinate of constituent i in the core relative 
to the core center-of-mass (see Fig|!J|). 

The M^EX) operator in Eq. (|Al[) can then be rewrit- 
ten as: 

.4-2 

M,(E\) = e^Zilrc + r'^Yx^rc + r'A 



2 



3=1 



(A2) 
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where the index i runs over the A — 2 constituents in the 
core, and j labels the two external nucleons. 

Making use of Eq.(l) in ref.[62| (see also [63|]) one can 
easily see that for A = 1 one has: 

3 

M^El) = eY,Z l r l Y 1 .^r l ) + M^El,core) (A3) 



APPENDIX A: TRANSITION STRENGTH 
OPERATORS 

For a system made of A constituents, the electric tran- 
sition operator of order A, M l _ l (EX), is defined as: 



A 



M„(EX) = e^Zi^Yx^h) 



where now i runs over the three constituents, Z; refers to 
the charge (in units of e) of each of the three constituents, 
and M^(El,core) — e^- =1 Zir' i Yi :ll (f' i ) is the electric 
dipole transition operator of the core. 
In the same way, for A = 2 one finds: 



(Al) M fl (E2) = eJ2z i r^Y 2 ^(f i )+M^E2,core) (A4) 
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^ f(m, n)r c Yi tm (r c )Mn- m (El,core) 



whose two first terms are analogous to the ones in 



Eq. (|A3p . The last term represents a coupling between 
the intrinsic core dipole operator and the dipole oper- 
ator associated to the center-of-mass of the core while 
f(m, /i) is a well defined function of the indices m and /i. 
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